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Abstract — The  performance  of  a  parametric  space-time 
adaptive  processing  (STAP)  method  is  presented  here. 
Specifically,  we  consider  signal  detection  in  additive  dis¬ 
turbance  containing  compound-Gaussian  clutter  plus  ad¬ 
ditive  Gaussian  thermal  white  noise.  Performance  is 
compared  to  the  normalized  adaptive  matched  filter  and 
the  Kelly  GLRT  receiver  using  simulated  and  measured 
data.  We  focus  on  the  issues  of  detection  and  false  alarm 
probabilities,  constant  false  alarm  rate  (CFAR),  robust¬ 
ness  with  respect  to  clutter  texture  power  variations, 
and  reduced  training  data  support. 

I.  Introduction 

This  paper  undertakes  a  performance  comparison 
of  several  candidate  space-time  adaptive  processing 
(STAP)  algorithms  in  compound-Gaussian  clutter  for 
airborne  radar  applications.  The  STAP  problem  is 
equivalent  to  hypothesis  testing  on  a  complex  (base¬ 
band)  measurement  (test  data)  vector  x  6  CJNwith  J 
channels  and  N  time  pulses.  Typically,  x  contains  an 
unwanted  additive  disturbance  d  with  unknown  covari¬ 
ance  matrix  and  may  contain  a  desired  signal  ae 
with  unknown  complex  amplitude,  a,  and  known  signal 
steering  vector  e.  The  binary  detection  problem  is  to 
select  between  hypothesis  Hq  :  a  =  0  and  Hi  :  a  ^  0 
given  a  single  realization  of  x. 

Current  research  [1-10]  is  addressing  the  detection 
problem  wherein  d  contains  partially  correlated  clutter 
described  by  a  product  model  [11].  Here,  the  clutter 
is  modeled  as  a  Gaussian  process  with  random  power 
variations  (scale  changes)  over  range.  This  model  is 
the  basis  of  the  spherically  invariant  random  process 
(SIRP)(or  compound-Gaussian)  clutter  model,  which 
includes  the  Weibull  and  K-distributions  as  special 
cases. 

In  [6, 12, 13]  the  optimal  processor  for  detecting  a 
rank  one  signal  in  SIRP  clutter  was  shown  to  be  equiv¬ 
alent  to  a  matched  filter  compared  to  a  data  dependent 
threshold.  With  a  simple  normalization,  this  test  can 
be  cast  in  the  form  of  the  normalized  matched  filter 
(NMF)  test  compared  to  a  data  dependent  threshold, 
the  calculation  of  which  requires  knowledge  of  the  un¬ 
derlying  clutter  probability  density  function  (PDF).  De¬ 
termination  of  the  clutter  PDF  imposes  onerous  train¬ 
ing  data  requirements.  Consequently,  ad  hoc  methods 
have  been  popular  in  recent  studies  [7-10].  A  popular 


ad-hoc  method  is  the  NMF  test  compared  to  a  data  in¬ 
dependent  threshold,  which  was  independently  derived 
in  [7,8]. 

The  work  of  [14]  (and  references  therein)  considered 
an  invariance  framework  for  hypothesis  testing  in  Gaus¬ 
sian  noise  having  a  covariance  matrix  with  known  struc¬ 
ture  but  unknown  level.  Interestingly,  the  test  statistic 
reported  in  [14]  is  identical  to  the  NMF  of  [7,8].  The 
work  of  [14]  also  extended  the  NMF  test  to  allow  for  an 
unknown  noise  covariance  matrix  and  unknown  scal¬ 
ing,  r?2,  denoting  the  ratio  of  the  test  and  training  data 
variances.  We  refer  to  this  test  as  the  normalized  adap¬ 
tive  matched  filter  (NAMF).  The  invariance  principle  of 
[14]  (and  references  therein),  and  perforce  constant  false 
alarm  rate  (CFAR),  applies  only  when  all  the  training 
data  vectors  are  scaled  identically  [10]. 

In  SIRP  clutter,  however,  each  training  data  vector 
realization  is  scaled  by  a  different  random  parameter. 
Although  the  NAMF  has  no  known  optimality  prop¬ 
erties  for  SIRP  clutter,  it  has  the  important  feature 
of  minimizing  dependence  upon  texture  power.  Some 
performance  results  of  the  NAMF  in  SIRP  clutter  are 
presented  in  [7,8]. 

Multichannel  model-based  (i.e.,  parametric)  methods 
for  target  detection  and  estimation  in  clutter  are  de¬ 
scribed  in  [2-4,  9, 10, 15-17].  In  particular,  a  model- 
based  STAP  method  called  the  non-Gaussian  para¬ 
metric  adaptive  matched  filter  (NG-PAMF),  requiring 
knowledge  of  the  underlying  clutter  statistics  was  first 
proposed  in  [3]. 

In  this  paper,  the  performance  of  the  normalized 
parametric  adaptive  matched  filter  (N-PAMF)  [10, 18]  is 
evaluated  and  compared  with  that  of  several  candidate 
STAP  algorithms.  Its  form  is  the  model-based  approx¬ 
imation  of  the  NAMF.  Statistical  equivalence  of  the  N- 
PAMF  test  to  the  NG-PAMF  test  was  shown  in  [10]. 
However,  unlike  the  NG-PAMF,  the  N-PAMF  test  re¬ 
quires  no  'a  priori’  knowledge  of  the  disturbance  statis¬ 
tics  [10].  This  feature  is  important  in  real-time  applica¬ 
tions  where  such  information  is  lacking.  Robustness  of 
Pd  over  a  broad  range  of  K-distribution  shape  parame¬ 
ters  (a)  ranging  from  Gaussian  (a  =  oo)  to  high-tailed 
PDF  (a  =  0.1)  is  presented  here.  These  considerations 
enable  assessments  of  CFAR  performance  with  respect 
to  the  amplitude  probability  density  function  (APDF) 
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associated  with  clutter  texture  variations.  Finally,  we 
examine  performance  versus  data  support  size  used  for 
disturbance  estimation.  This  issue  is  of  considerable 
importance  in  applications  where  training  data  support 
is  limited. 


II.  The  Clutter  Model 


A.  Non- Adaptive  Test  Statistics 

For  known  R^,  the  optimal  test  for  detecting  a  rank 
one  signal  in  Gaussian  interference  is  given  by 


jefRjV 

eHRrf 1 


H 1 

<  A mf- 

H0 


(4) 


Clutter  observed  in  a  single  channel  admits  a  repre¬ 
sentation  of  the  form 

Cfc(n)  =  vk(n)gk(n)  (1) 

where  a  complex-Gaussian  process  gk(n)  (speckle  com¬ 
ponent)  corresponding  to  time  n  and  range  cell  k  is 
modulated  by  a  statistically  independent  non-negative 
process  vk(n)  (texture  component).  Frequently,  ujb(n) 
is  approximated  as  a  random  variable  V  over  &,  but 
constant  over  time  if  it  has  long  temporal  coherence. 
Thus,  (1)  reduces  to  the  representation  theorem  [11] 
for  an  SIRP  and  is  given  by 

cjfe(n)  =  vkgk(n).  (2) 

For  the  multichannel  problem,  the  clutter  in  each  of  the 
J  channels  is  given  by  (2).  The  PDF  of  V,  fv(v),  is  de¬ 
fined  to  be  the  characteristic  PDF  of  the  SIRP.  The 
amplitude  of  ck(n)  is  K-distributed  for  Generalized- 
Chi  distributed  fv(v)  and  includes  the  Gaussian  model 
(a  =  oo)  as  a  special  case.  The  disturbance  d  con¬ 
tains  partially  correlated  clutter  c  modeled  by  a  Re¬ 
distributed  amplitude,  with  PDF 

OOL- f  lrQ 

Mr)  =  JL—-^Ke,-l(0r)  r  >  0,  0,  a  >  0  (3) 

where  0  and  a  are  the  distribution  scale  and  shape 
parameters,  respectively,  Ku(.)  is  the  modified  Bessel 
function  of  the  second  kind  of  order  i/,  and  T(.)  is  the 
Eulero- Gamma  function.  Small  values  of  a  result  in 
heavy-tails  for  the  PDF  of  (3).  From  (2),  the  clut¬ 
ter  covariance  matrix  is  Rc  =  R gE(V2)  where  Rg  6 
C JNxJN  jg  COvaxiance  of  the  Gaussian  (speckle) 

component  and  E(V2)  relates  to  the  texture  power. 

In  practice,  R d  is  unknown,  and  must  be  estimated 
from  a  signal-free  JN  x  K  secondary  data  matrix,  Z, 
whose  columns  are  assumed  to  be  statistically  indepen¬ 
dent  and  identically  distributed  (iid)  as  the  test  data. 
For  Gaussian  disturbance,  the  maximum  likelihood 
(ML)  estimator  is  the  sample  matrix  Rd  =  ZZH /K. 
However,  Rd  is  not  the  ML  estimate  for  compound- 
Gaussian  clutter. 


In  some  instances,  the  test  data  vector  can  have  a  co- 
variance  matrix  7?2Rd,  where  rf  is  an  unknown  level. 
The  phase  invariant  matched  filter  (PI-MF)  test  for 
these  problems  is  expressed  as  [14] 


A PXMF  ~ 


|egR/xj2 

T)2eH  RJxe 


Hi 

<  A PIMF 
H0 


(5) 


where  e  and  x  are  the  concatenated  JN  x  1  signal 
'search’  steering  and  data  vectors,  respectively.  The 
inner  product  of  whitened  vectors  b  =  Rd  2x  and 


f  =  RJ¥  e  is  the  matched  filtering  operation.  Although 
(5)  does  not  require  knowledge  of  signal  phase,  it  does 
require  knowledge  of  the  level  rj2  to  be  CFAR.  The  op¬ 
timal  test  for  this  problem  is  the  NMF  test  [14]  given 
by 


A  nmf  = 


|e*RdLxl2 

[e*Rjle][xHR^x] 


Hi. 

<  A  NMF- 
H0 


(6) 


The  test  statistic  of  (6)  is  simply  the  squared  magnitude 
of  the  inner  product  of  the  vectors  f  and  b  normalized 
by  their  squared  norms,  so  that  0  <  A  nmf  <  1. 


B.  Adaptive  Test  Statistics 

For  the  adaptive  problem,  R<f  replaces  Rd-  Conse¬ 
quently,  the  adaptive  version  of  the  test  of  (4)  denoted 
as  the  AMF  is  given  by 


|e"RIx*|2 
e^RJ1 


Hi 

<  A amf- 

H0 


(7) 


Observe  that  A  amf  is  simply  the  adaptive  version  of 
Amf -  Fo r  the  special  case  of  7}  =  1,  this  test  was  devel¬ 
oped  independently  in  [19,20]  where  its  CFAR  behavior 
was  noted.  This  property  is  lost  when  r)  ^  1. 

The  adaptive  version  of  the  test  of  (6)  is  given  by 


An  amf  = 


|egRJ1x|2 


[e"Rd  xx]  «o 


<  A namf ■  (8) 


Another  adaptive  detection  test  known  as  the  Kelly 
GLRT  [21]  is  expressed  as 


III.  Test  Statistic  Descriptions 

We  now  consider  several  non-adaptive  and  adaptive 
detection  test  statistics  in  this  section. 


A glrt  — 


|egRJ1x|2 


[eHR-le][l+_^_] 


<  KXglrt  (9) 

H0 


K 
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where  0  <  A glrt  <  1.  For  K  — >  oo,  the  tests  of  (7) 
and  (9)  converge  to  the  test  of  (4),  whereas  the  test  of 
(8)  converges  to  that  of  (6). 

In  this  paper,  we  consider  the  performance  of  the 
tests  of  (7)-  (9)  in  compound- Gaussian  clutter.  No  op¬ 
timality  or  CFAR  claims  of  these  tests  can  be  made  for 
the  case  of  SIRP  disturbance. 


C.  Model-Based  STAP  Tests 

For  multichannel  model-based  methods  [15],  the 
whitening  operation  is  performed  using  prediction  error 
filters  (PEF)  involving  time  series  or  state  space  archi¬ 
tectures.  We  define  y p(n)  as  the  Jxl  vector  error  resid¬ 
ual  output  of  a  P^-order  multichannel  linear  filter.  For 
a  multi-channel  (vector)  autoregressive  model,  a  tapped 
delay  line  architecture  is  used  where  the  Pth  order  filter 
coefficients,  A(&),  are  estimated  from  Z  using  a  multi¬ 
channel  parameter  estimation  algorithm.  These  J  x  J 
matrix  coefficients  provide  a  succinct  description  of  the 
spatio-temporal  disturbance  correlation.  Specifically, 

yp(n)  =  D0‘5L“1up(n) 

p 

=  D0  4Lo  A(*)x(n-  k  +  P|tfi)] 

n  =  0, 1, .  J V  -P-1 

(10) 

where  (10)  implicitly  defines  the  temporally  whitened 
Jxl  error  residual  up(n),  with  covariance  In 
practice,  Eu  is  unknown  and  the  estimation  algorithms 
produce  its  estimate  Eu.  However,  this  paper  em¬ 
ploys  obtained  by  averaging  the  outer  products  of 
up(n),  where  Up(n)  results  from  a  transformation  of 
the  secondary  data  by  the  prediction  error  filter  with 
fixed  A (k).  The  LDLH  decomposition  of  Ew  yields 
(L0,D0)  wffiich  are  used  to  spatially  whiten  u p(n)  [15]. 
Finally,  yp(n)  is  obtained  by  a  transformation  of  x(n) 
by  the  multichannel  prediction  error  filter  with  coeffi¬ 
cients  (Lo,Do)  and  A (k)  as  denoted  in  (10).  Similarly, 
the  transformed  steering  vector  s(n)  is  obtained  by  se¬ 
quencing  the  sequential  form  of  the  ‘search5  steering 
vector  e(n)  through  the  PEF  [3, 16].  The  normalized 
parametric  adaptive  matched  filter  (N-PAMF)  [10, 18] 
is  now  defined  as 


A  N-PAMF 


N-P-l 

]C  sH(-n)yp(n) 

71—0 

2 

'N-P-l 

J2  sH(n)s(n) 

.  n=0 

E  y?Wy p(n) 

.  n=0 

(11) 

A  related  parametric  adaptive  matched  filter  (PAMF) 
test  was  first  derived  in  [3]  for  Gaussian  disturbance. 


The  PAMF  test  is  identical  to  (11)  but  excludes  the 
second  bracketed  denominator  term.  In  [16],  several 
estimation  algorithms  are  considered  in  the  PAMF  im¬ 
plementation  for  Gaussian  disturbance.  In  this  paper 
the  multichannel  least  squares  method  is  used  for  filter 
parameter  estimation. 


IV.  Analytical  Results 

The  probability  of  false  alarm  and  probability  of  de¬ 
tection  for  NMF  operating  in  compound-Gaussian  clut¬ 
ter  (without  background  white  noise)  are  given  by  [9] 

Pfa-NMF  =  P(AjVMF  >  ^NMf\Ho)  =  (I-Aatmf)^"1 

(12) 

Pd— NMF  =  1  —  (1  ~  A NMF)JN  1 

*fMT^)Vu-W)]  <13> 

T(JN) 

where  g(.)  is  defined  in  [9]  and  6*  =  r(fc  +  1)r(Jjv  _  fc)  - 

The  expressions  for  Pfa  and  Pd  for  the  NAMF  oper¬ 
ating  in  Gaussian  clutter  are  given  by  [9] 


Pfa- NAMF 


/r(7) 


[1  +  (1  —  7)A  namf}4 


‘d—NAMF 


x  [1  —  gammainc 


-f 

=  1-/‘ _ i _ 

Jo  [1  +  (1  -  7)A namf]l 

E  (  ^ 

m-1  X  7 

4(1-7) 


-  dr,  (14) 


(r 


+  (1  -  7)  A  NAMF 


,mj]/r(7)d7 

(15) 

where  T,  gammainc(.)  and  /r(7)  axe  defined  in  [9]. 

Corresponding  Pfa  and  Pd  expressions  for  the  NAMF 
operating  in  SIRP  clutter  are  difficult  to  derive.  This  is 
due  to  the  fact  that  the  ML  estimate  of  the  SIRP  covari¬ 
ance  matrix  is  not  available  in  closed  form  [22].  Hence, 
the  NAMF  PDF  cannot  be  determined  analytically. 

The  Pfa  and  Pd  expressions  for  the  N-PAMF  and 
PAMF  operating  in  both  Gaussian  and  non-Gaussian 
clutter  scenarios  are  lacking  since  the  analysis  be¬ 
comes  mathematically  intractable.  Consequently,  per¬ 
formance  evaluation  of  these  methods  is  carried  out  by 
Monte-Carlo  techniques. 


V.  Performance  Results 

Performance  is  now  presented  for  the  detectors  de¬ 
scribed  above.  Probability  of  detection  (Pd)  is  com¬ 
puted  for  a  Pfa  of  0.01  via  100,000  Monte- Carlo  tri¬ 
als  using  a  physical  model  of  an  airborne  radar  sce¬ 
nario  [4].  The  target  signal  is  located  at  a  normal¬ 
ized  Doppler  frequency  fdt  ~  0.15  (unless  otherwise 
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Fig.  1.  Pd  versus  output  SINR  in  Gaussian  disturbance 

stated)  and  azimuth  0  =  0.  The  clutter  ridge  is  located 
along  the  normalized  angle-Doppler  plane  diagonal  with 
a  40  dB  (per  pulse,  per  channel)  clutter-to-noise  ra- 
tio(CNR).  The  one-lag  clutter  temporal  correlation  pa¬ 
rameter  [15]  is  0.999.  Disturbance  correlation  estimates 
are  obtained  using  K  secondary  data  cells.  The  output 
signal-to-interference  plus  noise  ratio  (SINR)  is  defined 
as  SINR  =  |a|2e/fR^1e.  Ail  performance  results  are 
obtained  with  compound-Gaussian  clutter  plus  additive 
thermal  white  noise. 

Figure  1  shows  Pd  versus  output  SINR  for  Gaussian 
disturbance  with  rp  ~  1.  Shown  here  are  analytical 
Pd  plots  for  the  MF  and  NMF  with  known  Rd.  The 
analytical  NAMF  Pd  curve  is  also  shown  for  K=128. 
Monte-Carlo  generated  Pd  results  for  the  PAMF(MLS), 
N-PAMF(MLS),  NAMF,  and  Kelly  GLRT  are  depicted. 
Performance  of  the  Kelly  GLRT  and  the  NAMF  are 
indistinguishable  for  this  example.  Note  that  the  N- 
PAMF  method  with  P  =  3,  nearly  achieves  the  NMF 
performance  for  low  sample  support  size  K=12.  Sin¬ 
gularity  of  Rd  for  K=12  precludes  implementation  of 
the  AMF,  NAMF,  and  Kelly  GLRT.  Figures  2  and  3 
display  Pd  versus  output  SINR  for  the  NMF,  NAMF, 
N-PAMF(MLS),  Kelly  GLRT  and  AMF  receivers  for 
clutter  processes  with  shape  parameters  a  =  0.5  and 
a  =  0.1,  respectively.  Observe  the  robust  behavior 
of  the  N-PAMF,  NAMF,  Kelly  GLRT  in  compound- 
Gaussian  clutter.  The  Kelly  GLRT  outperforms  the 
NAMF  at  high  SINRs,  whereas  this  condition  reverses 
at  low  SINRs.  Figure  4  plots  Pd  versus  the  clut¬ 
ter  shape  parameter  a  at  output  SINR=6dB  with  a 
ranging  from  0.1  to  1,000.  For  the  K-distribution, 
a  >  4  approximates  the  Gaussian  case.  The  results 
reveal  the  robustness  of  the  N-PAMF  and  NAMF  tests 
over  a  wide  range  of  shape  parameters.  However,  the 
N-PAMF(MLS)  shows  superior  performance  approach¬ 
ing  that  of  the  NMF.  Performance  of  the  PAMF  and 
AMF  degrade  in  heavy-tailed  compound-Gaussian  clut- 


Fig.  2.  Pd  versus  SINR  in  K-distributed  Clutter  a  =  0.5 


Fig.  3.  Pd  versus  output  SINR  in  K-distributed  Clutter  a  —  0.1 

ter.  However,  for  a  >  100  (Gaussian  region),  they 
incur  no  performance  degradation.  Figure  5  shows 
Pfa  versus  shape  parameter  a  with  each  test  statistic 
threshold  held  fixed  to  obtain  Pfa  =  0.01  for  Gaus¬ 
sian  disturbance  (a  =  oo).  A  significant  increase  in 
Pfa  for  the  NAMF  and  Kelly  GLRT  confirms  their  lack 
of  CFAR  with  respect  to  texture  variations.  Figures 
6,  7,  and  8  depict  plots  of  P/a  versus  threshold  for 
the  Kelly  GLRT,  NAMF,  and  N-PAMF,  respectively, 
for  several  K-distribution  shape  parameter  values.  The 
curves  for  the  N-PAMF  exhibit  much  lower  variability 
compared  to  the  Kelly  GLRT  and  NAMF,  reflecting 
a  robust  CFAR  performance  with  respect  to  the  tex¬ 
ture  PDF.  Figures  9  and  10  plot  the  test  statistic  vs 
range  cell  using  data  from  the  Air  Force  Research  Labo¬ 
ratory  (AFRL)  Multichannel  Airborne  Radar  Measure¬ 
ment  (MCARM)  program  with  an  inserted  target  signal 
at  range  bin  index  310.  Specifically,  data  correspond¬ 
ing  to  acquisition  ‘220’  from  flight  ‘5’  cycle  ‘e’  is  used  in 
the  examples  presented  here.  For  these  results,  we  use 
J=8  and  N=32.  We  define  the  performance  measure  $1 
as  the  ratio  of  the  test  statistic  at  the  test  cell  to  the 
mean  of  the  test  statistics  formed  from  adjacent  cells, 
and  $2  as  the  ratio  of  the  test  statistic  at  the  test  cell 
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Fig.  4.  Pd  versus  shape  parameter  (a) 


Shape  Parameter(a) 

Fig.  5.  Pja  versus  clutter  shape  parameter(a):frxed  threshold 


to  the  highest  test  statistic  formed  from  adjacent  cells. 
Figure  9  plots  the  test  statistics  for  the  NAMF  with 
K=512  and  N-PAMF(MLS)  (P=2)  with  K=16.  Fig¬ 
ure  10  plots  the  test  statistics  for  the  Kelly  GLRT  with 
K=512  and  the  N-PAMF  (MLS)(P=2)  with  K=16. 
Table  1  shows  $1  and  $2  for  several  values  of  K  and  P 
using  the  N-PAMF.  Note  that  the  N-PAMF  with  P=2 
and  K=16  provides  the  best  performance  for  this  sce¬ 
nario.  In  this  instance,  large  sample  support  does  not 
result  in  improved  performance  due  to  training  data 
nonhomogeneity. 


P 

(dB) 

*2  (dB) 

NAMF  (K=512) 

16.2 

7.45 

Kelly  GLRT  (K=512) 

16.71 

8.03 

N-PAMF  (K=512) 

4 

19.3 

12.3 

N-PAMF  (K=32) 

3 

22.1 

14.7 

N-PAMF  (K=16) 

3 

21.6 

15.2 

N-PAMF  (K=32) 

2 

21.8 

11.3 

N-PAMF  (K=16) 

2 

22.4 

14.7 

Table  1:  Values  of  and  for  the  N-PAMF(MLS), 
NAMF,  and  Kelly  GLRT 


Fig.  6.  Pja  versus  threshold  (Agljjt)  lor  the  Kelly  GLRT 


Fig.  7.  Pja  versus  threshold  for  the  NAMF 


Fig.  8.  Pja  versus  threshold  for  the  N-PAMF 


Fig.  9.  Test  Statistic  versus  Range 
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Fig.  10.  Test  Statistic  versus  Range 


VI.  Summary  and  Future  Research 

In  this  paper,  we  have  evaluated  the  performance  of 
the  N-PAMF,  NAMF,  and  Kelly  GLRT  in  compound- 
Gaussian  clutter  disturbance.  The  robust  detection  per¬ 
formance  of  these  methods  was  shown  over  a  wide  range 
of  clutter  texture  power  variations  (shape  parameters) 
for  K-distributed  clutter  processes.  Performance  of  the 
N-PAMF  was  found  to  be  close  to  the  NMF.  Next,  the 
CFAR  behavior  was  considered  by  observing  the  detec¬ 
tion  threshold  variations  with  respect  to  shape  param¬ 
eter.  Additionally,  we  demonstrate  the  robustness  of 
the  N-PAMF  method  with  respect  to  small  sample  sup¬ 
port  size  K  (secondary  data  cells)  used  to  estimate  the 
disturbance  correlation.  This  property  is  significant  in 
operational  scenarios  involving  range  varying  nonsta¬ 
tionary  clutter  which  severely  limits  the  availability  of 
representative  training  data.  Examples  with  real  data 
illustrate  the  potential  for  considerable  performance  im¬ 
provement  of  the  N-PAMF  over  the  NAMF  and  Kelly 
GLRT. 
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